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ABSTRACT 

The development of a high average power FEL for military applications would 
represent a significant improvement in missile defense, especially shipboard self-defense. 
The LANL regenerative amplifier FEL (RAFEL) is designed to produce an average output 
power of 1 kW. This FEL represents a significant increase in average power 
demonstrated in an FEL and also provides a test of the concept of combining the FEL 
oscillator and amplifier designs. Simulations were performed to better understand the 
physics behind the LANL RAFEL operation. 

Simulations study the transverse effects due to optical guiding by the intense 
electron beam and feedback. These simulations are applied to optimizing the undulator 
taper rate, feedback optimization, and initial phase velocity. Additional simulations study 
the longitudinal effects due to short electron pulses and optical pulse development over 
multiple passes. Finally, simulations of the RAFEL design using an ideal beam expand on 
understanding of the design’s basic characteristics and limitations. 
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I. INTRODUCTION 



A. INADEQUATE SHIP SELF-DEFENSE 

The world we live in today is very complex and dangerous. At the beginning of 
the decade the United States had a single adversary and a well defined mission. Since the 
end of the “cold war” both this mission and this single adversary have gone away. On the 
surface this sounds like a good thing. In reality it is a very unstable situation. The tight 
control over nuclear weapons and missiles is no longer there. Consequently missile 
technology and nuclear weapons are available to any country or organization with enough 
money to buy them. Additionally, the focus of the United States Navy is no longer that of 
a “blue water” navy. We no longer plan to fight large naval battles in the middle of the 
ocean. The mission of late has been in the littoral environment. Operating our ships near 
land puts them in the line of fire of even the crudest missiles. The close proximity to fire 
and proliferation of missile and nuclear technology make a good case for enhancing ship 
self-defense. 

Possibly more significant is the rate at which offensive missile technology 
continues to advance. Current approaches at ship self-defense are unable to keep up. The 
Close-In Weapons System (CIWS) is ineffective. The system just can not hit the newer 
missiles at sufficiently large ranges to protect the ship from debris. If the CIWS hits the 
missile at close range, the fragments would continue onward to inflict serious damage on 
the ship. The defensive missile concept has serious drawbacks. In this era of significant 
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cutbacks the excessive cost of each missile is unacceptable. In order to effectively defeat 
newer offensive missiles, new defensive missile cost is skyrocketing. This brings us to the 
inherent flaw in the defensive missile concept. If technology is invented for the defensive 
missile to improve its design over the offensive missile that technology is available for the 
offensive missile. This is unacceptable because the defensive missile has to be significantly 
better than the offensive missile to defeat it. A new weapon is needed as an alternative; 
this weapon is a “speed of light” weapon. 

This type of weapon would have speed of light response and could engage a target 
on the horizon. This is significant considering the prevalence of nuclear weapons we are 
facing. Destroying one on the horizon is infinitely more desirable than a few hundred 
yards one might get with CIWS. The concept of using a directed energy weapon, or a 
laser, for defense is not new. In fact, “technologies already exist and deployable HEL 
(High Energy Laser) weapons systems are just beyond the horizon, whether the United 
States Navy is pushing or not.” [1] 



B. THE FREE ELECTRON LASER (FEL) SOLUTION 

There are many directed energy and laser technologies developed and being 
developed. In fact, the MIRACL laser, a chemical laser, has been tested in White Sands, 
NM to shoot down various missiles. A few significant obstacles are power consumption, 
availability of space, and propagation through the atmosphere. 
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The power consumption problem depends on the choice of high energy weapon. 
For instance, a chemical laser would have to carry significant amounts of chemical reactant 
while a free electron laser (FEL) would need about 10 MW of power to operate 
continuously [2]. Substantially less power is nedded for operation in a limited 
engagement. This leads to the next major shipboard problem of where to put such a 
system on a ship. There just is not enough room on current ships, such as the DDG 5 1 , to 
add another weapon system. However, ref. [2] describes how an FEL and its prime power 
could be retrofit on a DDG 51. It could fit into the space of a 32 cell Vertical Launch 
System (VLS) weapon. When looking to minimize space and power consumption it is 
important to look at an HEL systems efficiency. The FEL may be able to achieve 10% 
system efficiency. However, a crucial problem remains, atmospheric propagation. 

The MIRACL laser was tested over a desert area. There are only small 
wavelength windows where a high energy laser beam can be successfully propagated over 
long distances. While it is possible to find a laser that will work at these wavelengths, it is 
not possible to adjust these lasers to variable weather patterns around the globe. The FEL 
is tunable and designable in wavelength. Once it is designed for a certain wavelength, it 
can be adjusted to another wavelength for optimum propagation with relative ease. 
Primarily for this reason the FEL is ideal for development as a shipboard self-defense 
weapon system. 

This thesis discusses a relatively new design for the FEL. In the past, FEL designs 
were either oscillators, which relied on many passes to achieve gain, or amplifiers, which 
have very high gain in a single pass. Los Alamos National Laboratory (LANL) has built 
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an FEL which combines both these designs. The LANL design is called a regenerative 
amplifier FEL (RAFEL) and relies on only a few passes to achieve high power. 

After a review of basic FEL theory, a detailed description of the RAFEL will be 
presented. The contribution of this thesis is to present results of simulations of the 
RAFEL to better understand the physics behind its operation. By simulating the 
diffractive effects inside the undulator, the optimum taper rate and feedback where 
obtained. These were optimized to achieve maximum operating efficiency. Using the 
same simulation the initial phase velocity was varied to explore possible wavelength 
variations. Finally, the pulse evolution along the undulator was simulated in order to 
understand the effects of short pulses. A better understanding of the physics of the FEL 
interaction can lead to more compact and more efficient FEL systems in the future. 
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II. FREE ELECTRON LASER PHYSICS 



A. INTRODUCTION 

The basic mechanism behind the free electron laser (FEL) is that electrons will 
radiate energy coherently when bunched at optical wavelengths. The FEL makes use of 
this process by bending a high energy beam of electrons back and forth causing them to 
radiate in an undulator. 

The main components of an FEL are an accelerator, an undulator, and for the 
oscillator design an optical resonator. There are also focusing elements and the source of 
electrons, an electron gun. The basic configuration is shown in Figure 1. 



electron beam 




resonator 



Figure 1. Basic components of an FEL. 
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A beam of electrons is injected into the accelerator which speeds them to 
relativistic velocities. These relativistic electrons are focused through an undulator. The 
undulator consists of static magnetic fields alternating in direction, normally from fixed 
magnets, which bend the electron beam. The beam of electrons then radiates forward in 
the direction of motion of the electrons. In the oscillator configuration this optical power 
is stored in the resonator cavity and will build over several passes until saturation is 
reached. Saturation is marked by decreased gain and by the field amplitude reaching its 
maximum value. In the amplifier configuration, there is much higher gain and the optical 
field grows significantly in a single pass. The physics of the electron motion in the 
undulator is governed by classical electromagnetic theory as described in the next section. 

B. PENDULUM EQUATION 

To understand the forces in the undulator, the transfer of momentum and energy 
between a free electron and an electromagnetic wave must be evaluated. A helical 
magnetic field in the undulator is assumed for the following calculations. The forces 
acting on the electrons in the undulator are governed by the following equations: 




( 1 ) 




( 2 ) 




( 3 ) 
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The optical electric field is E s = T'fcosT,- sin 4^,0] and the optical magnetic field is 
B s = £[sin V P, cos T,0] , where ¥ = kz - to/ + <}> . The optical wavenumber is 



k = to / c = 2n/ X and <{> is the optical phase. The velocity of the electron is given by cp. 
The undulator field in a helical undulator is B = 2?[cos(fc 0 z),sin(fc 0 z),o] , where B is the 



magnitude of the undulator magnetic field and k 0 = 2n / X 0 is the undulator wave number 
[3]. Equations (1), (2), and (3) describe the conservation of momentum, energy transfer 
and define the Lorentz factor, respectively. Insert the optical and undulator fields into (1) 



pendulum equation, it is useful to understand the FEL resonance condition. 

Think of the electron and photon racing along the undulator. In Figure 2, the 
photon is represented by a grey wavelength and the black dot is the electron. One 
undulator wavelength is shown as a dark curve. Resonance occurs when one wavelength 
of light passes over an electron while the electron travels through one undulator 
wavelength. Recall that the electron is traveling at a speed p.c where c is the speed of 
light. The condition for resonance is given by 



and integrate to get P ± = -K(cos k 0 z, sin k 0 z,0) / y where K = eB / k 0 mc 2 is called the 
“undulator parameter”. Taking this result and inserting into (2) yields 



ymc 

where ^ = (& + k 0 )z -cat and y-dy/dt. In order to continue in the development of the 




( 4 ) 




( 5 ) 



where X is the light wavelength and X 0 is the undulator wavelength. 
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Figure 2. Electron-Photon race. 

Inserting p ± into (3) gives a relationship between y and (3 Z . This can then be 
combined with (5) to rewrite (4) in the familiar form of the pendulum equation, 

.. IKeKE 



C = ‘ 



y 2 m 



-cos! 






( 6 ) 



This equation describes the electron phase dynamics as they evolve through the undulator. 

The amount of time the average electron spends in the undulator is given by the 
length of the undulator divided by the average electron velocity. Since bunching of 
electrons occurs over a length significantly smaller than the undulator length this time is 
about the same for all electrons. In order to convert to dimensionless units, a 
dimensionless time x = fict / L & ct / L will be introduced, where L is the length of the 
undulator and pc is the average electron velocity. According to this definition 0 < x < 1 
along the entire undulator. Substitute dt = Ldx / c into (6) for the final form of the 
pendulum equation 



OO 

C = |cr| cos(<^ + <{)) , 



( 7 ) 
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where C, = d 2 C, / dx 2 and |a| = 4nNeKLE / y 2 me 2 is the magnitude of the dimensionless 



optical field. This equation governs the motion of the electrons in phase space. In order 
to arrive at this final solution it was assumed that the electrons are relativistic and that the 
change in the electron energy is relatively small along the undulator [3]. Having 
developed the equation describing the electron trajectory in the undulator the physics of 
the optical wave must be understood. 

C. WAVE EQUATION 

In the FEL, the light wave is amplified by stimulated emission. An equation 
describing the light wave dynamics can be developed starting from Maxwell’s wave 
equation. 



The vector potential is A(x,t) = E(x,t)(sm'¥,cos v V,0) / k and J ± (x,t) is the total beam 



and <j> are slowly varying in z and t, and complete the left side of (8). This allows the wave 
equation to be written as 




( 8 ) 



current. Define the complex electric field S{x,t) = E{x,i)e'^ x ''\ and rewrite the vector 
potential A(x, t) = S{x,t)ee' a Ik where a = kz - cor and e = (-/,0,0) . Assume that both E 




( 9 ) 



The current is the sum over all single particle currents, 




( 10 ) 
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Identify C,= k 0 z + a = (k 0 + k)z - to t and simplify the sum by using the weight factor 
p(x,/) = electron density at 3c . Then assuming the positive solution of L, the wave 
equation becomes 






l- 2 Ad Id 
U v + '*l 7z + ~cT,j) 



S{x,t) = -2nieKkp(x,t)( — ) 

\ y / (i,«) 



(id 



where ( ?) is the average over a sample including a number of electrons in one periodic 
section of phase space located at position (x, t). In (1 1), S{x, t ) is an envelope of the 
optical field driven by the current density p(x,/) . The term ( defines the degree of 
electron bunching in the volume (3c, t). 

Equation (11) must be made dimensionless to correspond to the development of 
(7). In addition to the dimensionless time x = tc / L , the coordinate transformation 
z = z-ct makes the frame of reference that of the electron beam. Applying the chain rule 
and after some algebraic manipulation [4], the wave equation can be written as 




where j = 8ti 2 Ne 2 K 2 1} p / y 3 mc 2 and a = |tf|exp(/<{>) . If the electron pulse is long, the 
spatial dependence can be dropped and a simplified wave equation, 

a = -;(«-*} , (13) 

is obtained. The solution for the linear undulator follows the same procedure except that 
K is replaced by K(j 0 {%) - J x {t)) , where £ = K 2 ll(\ + K 2 ) [3]. 
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D. (LOW GAIN) FEL OSCILLATOR DESIGN 



The pendulum equation and wave equations derived previously form a self- 
consistent Maxwell-Lorentz FEL theory. The dimensionless current density,/, couples the 
electron beam to the optical wave and is proportional to the change in the optical field. 

This theory is valid for strong or weak fields, a> n or a < 7t respectively, and high or low 
current, / > n or j< n respectively. [3] 

An FEL oscillator generally operates in the weak field, low current regime. This 
design makes use of multiple passes to achieve gain and saturation. One method of 
achieving multiple passes is by reflecting light inside a resonator cavity as Figure 1 shows. 
In this case, one of the mirrors would have to allow some fight out. Typically this is done 
with a hole in the center of the mirror or a partially transparent mirror. 

Low current density,/, means there is a small change in the field a over a single 
pass. In the limit of low current density the field a can be assumed to be roughly constant 
and the integral of the pendulum equation, (7), becomes simple and yields, 

v 2 =2|a|sin(; + <j)) + 2// 0 . (14) 

H 0 is a constant of motion resulting from the conservation of energy, and is called the 

Hamiltonian. The electron phase velocity is defined as v = C, = z|(fc + fc 0 )P z - &] • A weak 

field generally means that there is little change in the electron phase velocity over the 
undulator. Therefore, in weak fields with low current the solution to the pendulum 
equation, (7), can be given by [3] 

;(x) = <; 0 + v 0 x - %[cos(^; 0 + v 0 t) - cos(<; 0 ) + v 0 x sin(^ 0 )]+... (15a) 

v o 
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and 



v(t) = v 0 + — [sin(^ 0 +V 0 x)] + 



£o 

vi 



2 r i 

T — — (cos(2^ 0 + 2v 0 x) - cos(2<; 0 )) - 1 - v 0 x sin(^ 0 )cos(C 0 + v 0 x) 



, (15b) 



■f... 



where C, 0 is the initial electron phase, ao is the initial dimensionless optical amplitude and 
v 0 is the initial phase velocity. 

Near resonance small changes in the phase velocity are given by Av= AnNAy / y 
and changes in phase are given by AC , « kAz . These equations relate v and C, of (15) to 
the electron energy and the microscopic position of the electron in the beam. Phase space 
evolution is an ideal way to track electron dynamics and bunching. The electrons enter the 
phase space with (v 0 ,^ 0 ) and then evolve according to (14). The specific phase space 
path of an electron is given by the constant H 0 [3]. The initial phase velocity is 
determined from the resonance condition discussed earlier, where v 0 = 0 at resonance. 

Phase space is most simply understood by tracking the familiar motion of an 
undriven pendulum. The points (-7r/2,0) and (3rc/2,0) are unstable fixed extrema and 
correspond to a pendulum at the top of its arc. The point (n/2,0) is a stable fixed extrema 
and corresponds to a pendulum at the bottom of its arc. In other words if a pendulum was 
released with no drag and allowed to swing back and forth it would trace out a closed 
orbit. Suppose this pendulum could come infinitely close to vertical in its swing path back 
and forth. This would trace out a path called the separatrix. The separatrix separates 
open orbits from closed orbits and is indicated by the solid line in Figure 3. However, if 
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the pendulum was given enough initial force to propel it over the top in a circular orbit it 
would correspond to an open orbit, one outside the separatrix. Setting H 0 = 0 in (14) 

gives the equation of the separatrix. It is evident that the height of the separatrix is 4-y/joj . 

Phase space plots are particularly useful in determining FEL operation by following 
the evolution of the electrons. A single 27t section of the phase space is adequate because 
it is identical to others nearby when the coherence length of the optical wave is long [3]. 
When looking at a phase plot, it is important to realize that electrons near C, = n transfer 
energy to the optical wave and those near C,= 0 take energy away from the optical wave. 
The phase velocity for maximum gain is v 0 = 2.6 [3]. 



*** FEL Phase Space Evolution ** ,k 

i «1 a=3 . 14 v *2 . 6 

J o o 

0.2 



0 

0.1 

0 

-it/2 £ 3tt/2 Ox 1 



Figure 3. Electron phase space plot with low current and weak fields. 
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Figure 3 is a phase space plot of electrons through an undulator with a weak field 
and low current density created by solving the pendulum equation, (7), and the wave 
equation, (13). The electrons are injected into the undulator with uniformly spaced phases 
and initial phase velocity v 0 = 2.6. Normally there would be some spread in the phase 
velocities, but for this example it was assumed to be negligible. The initial position of the 
electrons is light grey and the positions get darker as the end of the undulator is reached. 
Electrons that initiate outside the separatrix stay in open orbits and those initially inside 
stay in closed orbits. The electrons bunch near C, = n, and therefore have significant gain, 
G, which is plotted to the right with the optical phase, <j> with respect to x . At the 
beginning of the undulator, there is little gain or phase change, but as bunching develops 
toward the end the gain and phase increase. 

For low current and weak fields, energy conservation can be used to calculate the 
gain by equating the radiation energy of a volume of light with the energy lost from a 
volume of electrons. Electrons uniformly distributed in phase that are monoenergetic have 
an average energy loss of ymc 2 ((v)~ v 0 ) / 4nN . The first-order terms from (15) average 
to zero, but the second-order terms allow the gain to be written as 

G(t) = jF 

where F = (r b I w 0 ) is the “filling factor”. Here r b is the electron beam radius and w 0 is 
the optical mode radius. The gain is defined as the fractional increase in power contained 
in the optical field. In general, the interaction within the FEL changes both the amplitude 



2-2 cos(v 0 x) - v 0 x sin( v 0 x ) 
v* 



(16) 
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and phase of the optical field. Using the expansions given by (15) the wave equation, 
(13), can be solved and in lowest order of ao and j given by 

2-2 cos( v 0 x) — v 0 x sin( v 0 t) 



and 





f ( 


c? 

II 

"■p 

J? 


1 + 7 




l V 




f 



2v£ 



)j 

)) 



(17a) 



<|>(t) = j 



2sin(v 0 x) - v 0 x(l + cos(v 0 x)) 



2v£ 



(17b) 



In this weak-field low current case, the gain is dependent only on the current density and 
the initial phase velocity [3]. 



E. (HIGH GAIN) FEL AMPLIFIER DESIGN 

The amplifier design usually makes use of high gain in a single pass. The 
components are similar to the oscillator design except that there is no resonator or 
feedback loop. In order to achieve high gain, amplifiers generally operate in a high current 
density regime with long undulators. 

As the current density is increased, j » 7t, the optical phase change and optical 
amplitude change are no longer small. The solution to the pendulum equation, (7), given 
by (15) is not valid in this regime. However, the pendulum and wave equations (7) and 
(13) are valid in the high gain regime so they can be solved numerically to explore the 
phase space evolution of the electrons. Figure 4 shows the evolution resulting from a 
current density of j = 100 and uniformly spaced electrons with an initial phase velocity of 
v 0 = 1.6, which corresponds to maximum gain in this case [3]. Like the low current 
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*** FEL Phase Space Evolution *** 

j =100 a =3.14 v =1.6 

J o o 




-ti/ 2 £ 3rt/2 0 t 1 



Figure 4. Electron phase space plot with high current and weak fields. 

density plot in Figure 3, there is little bunching near the beginning of the undulator, x = 0, 
and therefore little gain or phase change. In this example, the beam is near resonance so 
the optical phase begins to grow without the corresponding increase in optical amplitude. 
This causes the separatrix to shift toward C,= it and high gain begins. The high gain in this 
example causes the height of the separatrix to increase dramatically. Figure 4 shows the 
evolution of the separatrix from light grey to black as the electrons move along the 
undulator from x = 0 to x = 1 . 

In this high current, weak field regime, the wave and pendulum equations can also 
be solved analytically using perturbation theory. The expanded wave and pendulum 
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equations can be combined and integrated over all initial phases, 1 dC , 0 , to obtain the 
following equation governing the evolution of the optical field [3] 



a(x)= lj ^\ldx'x'F{x')e t-t') 



( 18 ) 



where F(x') = J dqf{q)e ,qz ' is the characteristic function of the distribution of the initial 
electron phase velocities,/^. In (18), reference to electron phases has been removed, but 
there is an average over the distribution of initial electron phase velocities. Equation (18) 
is valid in weak fields with high or low current density. The form of the solution can be 
found by taking successive derivatives. After finding the roots to the characteristic 
equation and applying initial conditions the solution, 




e (y/2) 1/3 (/ + V3)i/2 



+ e (j/2) m (i-Sh/2 



+ e 



Ail 2)' 



1/3 




(19) 



is found. This shows that for small x , there is little change in the optical field amplitude. 
This remains true during the bunching time, x < x B « (2 / y)' \ As soon as the electrons 
bunch there is exponential growth and high gain. In the high current limit, the fastest 
growing term from (19) dominates, and the solution to the wave equation, (13) is given by 



|fl(x)| « -fe 



a 0 A ;/2) I/3 v 3t/2 

3 

1/3 



( 20 ) 



<|>(x) « (j / 2) ,A x / 2 



For high current, the optical gain is given by 

G(x) = ie ( " 2|WA . (21) 

In order to obtain the solutions given in (18) and (19), perfect beam quality and initial 
phase velocity, v 0 = 0 , were assumed. 
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1. The Tapered Undulator 



As electrons give up energy to the optical wave and slow down they begin to drop 
out of the gain spectrum. This is a problem in the amplifier design since it is important to 
get the maximum energy out of the electrons in a single pass. It would be ideal if the 
electrons could be accelerated to stay in the gain spectrum. Resonance would then be 
maintained and the interaction between the electrons and the optical wave would be 
allowed to continue. 

There are two basic methods for replacing or compensating for the electron’s lost 
energy. They are evident if one looks at the resonance condition, (5). The electron losing 
energy means that y is decreasing. This means that either the “undulator parameter K, 
or the undulator wavelength, A, 0 , must decrease to maintain resonance. Decreasing the 
undulator wavelength along the z axis means decreasing the spacing between the magnets 
as the electron progresses from the beginning to the end of the undulator. This results in 
an artificial acceleration of 8 * -2nNAk 0 / , where N is the number of periods and X 0 is 

the undulator wavelength. Another method is to decrease K, which is directly 
proportional to the undulator field strength, along the z axis giving an artificial 
acceleration of 8 « -4nNK 2 AB / i?(l + AT 2 ) , where B is the magnetic field amplitude. 
Figure 5 demonstrates these concepts. The arrows represent magnetic field inside the 
undulator, where a shorter arrow means a smaller field. The untapered undulator is 
represented by A. B and C represent tapering the magnetic field and undulator 
wavelength respectively. With tapering included (7) and (13), the pendulum and wave 



18 




equations become 

OO O O / V 

C, = v = 5 + \a\ cos(^ + <J>) and a = - 7 (e"' ? ) , ( 22 ) 

where C, = I A; u (x')^t' + kz -(tit [3]. These equations are valid for weak and strong optical 

0 

fields with either high or low gain, but not when efficiency is large, r| « 10 % or more. 

Tapering can be advantageous when the artificial acceleration 8 exceeds the 
deceleration obtained without taper in strong optical fields. The desirable criteria for 
tapering is given by 

|a|>8>4|ar 2 >2n , (23) 

where the left side of this inequality is the condition for trapping electrons in tapered phase 
space. The middle describes the condition that the artificial acceleration exceeds the 
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deceleration. The right side is the condition that the taper operate in the strong field 
regime. 

Electrons in the tapered undulator evolve a bit differently in phase space than in the 
untapered undulator. Electrons with an initial phase of <^ 0 = 0 are accelerated by both the 
optical field, \a \ , and the artificial acceleration, 8 , and will eventually contribute little to 
the interaction between the optical wave and the electron beam. Electrons with an initial 
phase of <^ 0 = n see the forces of 8 and \a\ roughly cancel out leaving them trapped in 
closed orbits. From this perspective, it is easy to see that tapering works because 
electrons which give energy to the optical wave are trapped and those which take energy 
from the optical wave are accelerated away. 

Tapering the undulator leads to much higher efficiency and gain in strong optical 
fields, but has lower gain in weak optical fields. When the taper rate, or phase 
acceleration, 8 is too large for the FEL field strength, electrons will not remain trapped 
and efficiency is reduced. Figure 6 is a phase space plot of 1000 electrons in an initial 
strong field of a 0 = 40 that is amplified by a current j - 1. The initial phase velocity is 
zero and the taper rate is 8 = 107t which meets the criteria of the (23). From this figure it 
is evident that the untrapped electrons are spread fairly evenly in phase, and the trapped 
electrons are concentrated inside the separatrix in the area for gain. Electrons can be seen 
to leak out of the separatrix on the right side due to the growth in optical amplitude and a 
shifting phase, <j> . The oscillations in the gain and phase plots on the right are due to the 
revolution of the electrons inside the separatrix. Examination of Figure 6 reveals that 
about 50% of the electrons are trapped. As it turns out the efficiency of this tapered 
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*** FEL Phase Space Evolution *** 
j =1 a Q =40 v Q =0 5=1 On 




-n/2 C 3n/2 0 t 1 



Figure 6. Tapered undulator phase space plot with low current. 



undulator is about 3 times the natural efficiency without taper [3]. Increasing the current 
traps more electrons and the gain increases even more significantly. 

2. The Regenerative Amplifier FEL (RAFEL) 

LANL has demonstrated a combination of the oscillator and amplifier designs 
called the RAFEL shown in Figure 7. The device is designed to operate in the high 
current, high gain regime, but uses a small amount of feedback. The photoinjector injects 
the electrons into the linear accelerator and then they pass through a hole in the first 
mirror. The electrons in the undulator give up energy to the optical wave. The beam is 
then guided through a hole in the second mirror at the end of the undulator and is directed 
into a beam dump. Most of the optical wave passes through the same hole in the mirror. 
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but a small fraction is recirculated through a feedback loop to the beginning of the 
undulator to be amplified by the subsequent incoming electrons. While feedback is a 
characteristic of oscillator designs, the concept behind the RAFEL’s feedback is 
fundamentally different. The purpose is only to seed the FEL with a small amount of 
power so that the process does not start from noise. The high single-pass gain overcomes 
the small amount of feedback. This allows the RAFEL to reach saturation in only a few 
passes. Details of the RAFEL design will be discussed in the next chapter. 



Regenerative Amplifier FEL 



Feedback Loop 




Electron Dump 



Figure 7. Diagram of the Regenerative Amplifier FEL. 



F. 3D DIFFRACTION 

Any freely-propagating, coherent optical wave will spread due to natural 
diffraction as it passes through the undulator. The Rayleigh range, Jtw 0 2 / X, is defined as 
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the distance over which light of wavelength, X, initially doubles in transverse area. The 
initial area of the optical wave is given by tva>1 . It is important to ensure that the Rayleigh 
range is comparable to the length of the undulator, L, or FEL coupling could be reduced 
by diffraction of the wave away from the electron beam. This occurs if the optical mode 
becomes significantly larger than the electron beam. If the current density is high enough, 
this restriction on the Rayleigh range can be modified because the radiation wavefront is 
distorted, or guided, by the electron beam. 

FEL optical guiding can occur when the high-current-density electron beam 
continually focuses the optical light wave back into the electron beam. The optical phase 

evolution along the undulator is dfyj =(jr e 2 / 2 vVq )' 3 dz / 2L. This includes the “filling 
factor” discussed previously. Natural diffraction produces a phase shift d§ D = -dzX / 7 m ’ 2 
which is of the opposite sign. In order for the FEL to have optical guiding d§ } > d$ D . 

From this inequality, it is evident that either a high j or a large electron beam are necessary 
for effective guiding of the optical mode. [6] 

The FEL theory developed earlier can be extended to include transverse modes. 
The complex optical field envelope is still taken to be slowly varying along z. The 
derivation is essentially the same as before. The transverse coordinates can be normalized 
to the characteristic mode size, x^n / LX — » x and yyfnTTX — » y, where L is the length 
of the undulator and A is the optical wavelength. The longitudinal coordinate can be 
normalized to the slippage distance, {z-ct / NX) — > z, where c is the speed of light and N 
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is the number of periods in the undulator. Then using dimensionless time, x = ct / L, the 
parabolic wave equation can be written as 



4 



dr. 



-tV? + 4~ 



(x>yw 



(24) 



where V* = 8^ + 8^ [3]. This equation governs the dynamics of the optical wave in the z 
direction as before, and the V_j; operator describes diffraction in the transverse (x, >>) 



directions. The solution to (24) is given by 

a(x, y,z,x + Ax) = exp^- V* jo(x, y, z, x) - A x(yV c ) w _ t x +... 



(25) 



to order jAx 2 [3]. The evolution of the electron current is governed by the Lorentz force 
equation for each electron in the beam. The pendulum equation can be generalized to the 
form 



v(x, y, z - x, x) = £(x, y, z — x, x) = |a(x, y, z, x)| cos[^(x, y, z - x, x) + <)>(x, y, z, x)] (26) 

to complete the theory in 4 dimensions [3]. 



G. TRAPPED PARTICLE INSTABILITY 

When saturation is reached in a strong field, the height of the separatrix is large 
and the electrons tend to travel in closed phase space paths causing increased efficiency. 
While in these closed paths, some electrons become trapped in deep potential wells. This 
results in decreased bunching and decreased gain. This effect is termed the “trapped- 
particle instability.” It occurs when an FEL evolves to power levels beyond saturation. 
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The oscillations of the trapped electrons starts the instability and for electrons near 
phase C, 0 the motion is described by 

C(x) = C o +— sinfv/c) , (27) 

for the initial position (<^ 0 ,v 0 ) [3], The trapped-particle oscillation frequency, or 

“synchrotron” frequency, is given by v f = (|<r| 2 -8 2 ) . The taper parameter, 8 , is zero 
when the undulator is not tapered. Sidebands appear at v 0 ± v s away from the fundamental 
wavelength by AX / X = v s / 2nN . The shift is simply the ratio of the number of 
synchrotron oscillations to the number of periods in the undulator. 

H. SHORT PULSE EFFECTS 

When the electron beam is made up of short pulses, comparable in length to the 
slippage distance NX, the interaction between the optical wave and the electron beam will 
be dominated by short pulse effects. A series of short electron pulses generates a series of 
short optical pulses. In the oscillator design these optical pulses bounce between the 
resonator mirrors separated by a distance S. The reflected optical pulse arrives at the 
beginning of the undulator in time intervals, 2 S / c, and the entering electron pulses must 
be synchronized to arrive at the same times. Desynchronism, d, is the displacement 
between the electron and optical pulses at the beginning of the undulator. Desynchronism, 
d, is normalized to the slippage distance so that when d = 0 there is exact synchronism. 

The light pulse drifts away from the electron pulse over many passes when there is 
exact synchronism. This effect is due to the light pulse being distorted on each pass 
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because gain is preferentially added to the trailing edge of the pulse. This causes the 
effective speed of the pulse to slow and fall away from the bunched electrons. The end 
result over many passes is no gain. The desynchronization can be controlled by adjusting 
the path length or reducing the distance between mirrors. 

At small desynchronism, the optical power can be large enough to cause trapped- 
particle instability, creating a broad optical spectrum and a broad electron spectrum. The 
trapped-particle instability is identified by sharp spikes in the optical pulse. Since a small 
change in d causes a large change in steady-state power, the FEL is generally unstable at 
small values of desynchronism. The FEL becomes more stable for larger values of d in the 
operating range. The reduced coupling at larger d causes lower steady-state power, and 
therefore trapped-particle instability generally does not occur. The power spectrum is 
narrow because the optical pulse is smooth and longer than the electron pulse. [3] 

Another effect of short pulses is “limit cycle behavior” in the optical pulse. This 
behavior is represented by the power and optical pulse oscillating periodically over many 
passes. Limit cycle behavior occurs at values of d where a stable laser pulse configuration 
does not occur. The structure initiates on the trailing edge of the optical pulse and the 
desynchronism mechanism forces it through the pulse envelope. [3] 
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III. RAFEL DESIGN 



A. INTRODUCTION 

As mentioned previously, the LANL RAFEL is a combination of the oscillator and 
amplifier designs. The term regenerative amplifier was coined for this design because it is 
similar to conventional laser regenerative amplifiers. The optical feedback is only on the 
order of 10" 2 to 10' 5 . As a consequence of the electron beam pulse structure, shown in 
Figure 8, the optical feedback path must be synchronized with the incoming electron beam 
to ensure there is overlap. The advantage of this type of FEL design is that the mirrors 
never get exposed to high peak or average power. [6] 




Figure 8. Pulse structure of LANL RAFEL [6]. 
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B. LANL RAFEL SYSTEM CHARACTERISTICS 



1. The Electron Beam 

The rf power source provides an output greater than 20 MW peak for 30 
microseconds. The average power of the 1300 MHz klystron transmitter is 50 kW in 
bursts [6]. This powers a photoinjector/linac which creates a beam of electrons with an 
energy of 17 MeV, 300 A micropulses at approximately 100 MHz. The current 
macropulse repetition rate is only 1 Hz, but 60 Hz will be available when diode-array- 
pumped amplifiers are installed. Details of the pulse structure are shown in Figure 8. The 
duty factors for the micropulse and macropulse are on the order of 0.2% [6]. This 
combined with the 300 A peak current and 17 MeV electron beam leads to the predicted 
kW output power. Multiplying the beam energy by the true average current gives; 

17 x 10 6 F x 300 A x 0.002 2 « 20 kW . This value multiplied by the efficiency gives the 
output power. If an efficiency of 5% can be achieved, this would give an output of 1 kW. 
The relatively high current density could cause guiding of the optical mode which would 
be favorable to FEL coupling. However, the small electron beam has a small radius of 
only 0.18 mm could reduce this effect. 

2. The Undulator 

The undulator is L = 2 meters long with a X. 0 = 2 cm period. Each period contains 
four permanent magnets in a modified Halbach configuration [6]. The first meter of the 
undulator is untapered while the last meter has about a 30 % taper in the undulator field. 
The undulator gap in the first meter is 0.59 cm. In order to accommodate the growing 
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optical mode the taper is accomplished by increasing the gap between the magnets to 
effectively reduce the magnetic field that the electrons experience. This is accomplished 
by increasing the undulator gap to 0.95 cm by the end of the tapered section of the 
undulator. The undulator period remains constant throughout the undulator. Two plane 
focusing is accomplished with a rectangular groove cut on the pole face of the magnets 
and is equal strength in both directions. 

3. Feedback Loop 

The feedback loop consists of mirrors as shown in Figure 9. The annular mirrors 
are relatively flat with the upstream mirror having a 0.25 cm hole in the center and the 
downstream mirror having a 0.70 cm hole in the center. These holes allow the electron 
beam to pass through unimpeded. The downstream mirror has a larger hole to allow the 
portion of the optical mode centered near the electron beam to pass through. The outer 
portion of the mode is reflected back through spherical and cylindrical mirrors to seed the 




Figure 9. Feedback loop for LANL RAFEL [6]. 
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undulator. Therefore, in order for the design to function as anticipated, a small amount of 
the optical field in a donut shape feeds back enough fight to continually seed the undulator. 
The large-signal optical field concentrated around the electron beam provides outcoupling. 
The detuning length for the RAFEL design is longer than conventional FEL designs and 
must be within 6 mm of the 20 ps electron pulse length to observe lasing. 



C. SYSTEM PARAMETERS 

Table 1 is a summary of some significant electron beam and undulator parameters 
used for the simulations in the following chapter. The experiment was set up for a peak 
current of 300 A. A peak current of 370 A is possible in the future so some simulations 
were conducted at that current for comparison. The undulator parameter fisted in the 
table is given by K = eB rms X 0 / 2nmc 2 , where X a is the undulator wavelength or period. 
From the resonance condition presented earlier, X = X Q (l + K 2 ) / 2y 2 , the expected optical 
wavelength of the RAFEL is about 16 pm. 



Electron Beam Parameters 


Undulator Parameters 


Electron kinetic energy 


17 MeV 


Length 


2 m 


Micropulse length 


18 ps 


Period (wavelength) 


2 cm 


Beam radius 


0.18 mm 


Magnetic field amplitude 


0.7 T | 


Emittance 


5 mm-mrad 


Undulator parameter 


0.92 


Peak current 


300 - 370 A 


Magnitude of taper 


30.25% 


Fractional energy spread 


0.1% 


Taper starts at 


1 m 



Table 1. LANL RAFEL parameters. 
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IV. SIMULATIONS 



A. INTRODUCTION 

This chapter discusses simulations of the RAFEL. The parameters of the LANL 
experiment are explored in order to develop a better understanding of the physics behind 
the RAFEL design. The goal in exploring the LANL parameters is to examine the 
feasibility of reaching a kW and to optimize a few of the parameters. In support of this, 
the taper rate and feedback amount were optimized using a simulation including diffraction 
in three dimensions (x,y,t), but not taking into account longitudinal pulse effects and 
slippage. The same simulation was used to explore the maximum efficiency versus the 
initial phase velocity, v 0 . In order to enhance understanding of the RAFEL concept, 
simulations including pulse effects and multiple passes were performed without diffraction. 
Further simulations were performed without diffraction and using a perfect beam to 
idealize the full potential of the RAFEL concept. 

B. TAPER OPTIMIZATION 

Three dimensional (3D) simulations were used to explore the effects of diffraction 
in the undulator with varying taper rate 5 . The LANL RAFEL is constructed with the 
taper starting halfway through the undulator. Figure 10 shows the output of a 3D 
diffraction simulation with a peak current of 300 A and feedback of 10* 4 . The taper rate 
in this case is 5 = 12871 corresponding to the design taper listed in Table 1. Recall that the 
undulator length indicated by x goes from zero to one. 
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Figure 10. 3D RAFEL simulation with 8 = 128rc . 

The evolution of the optical mode, |a(x,x)|, shown in the upper-left window of 
Figure 10, is dominated by diffraction after the taper is turned on at x = 0.5. The intense 
electron beam provides guiding over the first meter of the undulator, x < 0.5, but 

O 

diffraction overcomes this focusing over the last meter of the undulator where <j)(x) < 0. 
The window in the top-center shows a cross-sectional view of the optical mode, |a(x, y)| , 
at the end of the undulator. The top-right window lists the dimensionless parameters for 
this simulation; a 0 = 3.3 is the initial dimensionless optical amplitude, j = 5.5 x 10 4 is the 
dimensionless beam current density, g x = cr v = 0.053 corresponds to the dimensionless 
radial size of the electron beam, and <r 6> , = 0.9 and ct c = 1.3 characterize beam quality. 
The random Gaussian spread of the electrons is a G and the spread in angle is a e>1 . The 
graph at the right-middle shows the development of the optical phase along the undulator 
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length <|)(t) . When the taper turns on at x = 0.5, the optical phase change along the 

O 

undulator length becomes negative 4>(x) < 0 and encourages the optical wavefront to be 
excluded from the electron beam, thus decreasing the interaction strength. The lower- 
right is the plot of the dimensionless optical power growth along the undulator P(x) and 
the plot of the natural logarithm of the single-pass gain along the undulator ln( 1 + G(x)) . 
The middle-left window shows the bunching current development in the electron beam 
along the undulator, with the end view in the center. This is indicated by the appearance 
and darkening of a line, indicating bunching current, in the center of the plot. At the 
bottom- left is the beam’s electron phase velocity evolution along the undulator, and the 
bottom-center is the beam’s electron phase space plot at the end of the undulator. 

In order to optimize gain and efficiency for an FEL, about 50% of the electrons 
should be trapped [4]. In Figure 10, much less than half the electrons are trapped near 
resonance v = 0 as can seen in the bottom-center photo where there is no significant 



bunching, indicating the taper rate is too large causing electron detrapping. 
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P(z) 


G f 


8 O 71 


7.1% 


2.4x1 0 5 


1.67xl0 3 


90ti 


7.7% 


2.5xl0 5 


1.77xl0 3 


95tt 


7.9% 


2 . 5 x 1 0 5 


1.80xl0 3 


10071 


8 . 1 % 


2 . 6 x 1 0 5 


1.83xl0 3 


1 1 Oti 


8.3% 


2 . 6 xl 0 5 


1.87xl0 3 


11271 


8.3% 


2.6xl0 5 


1.87xl0 3 


11571 


8.3% 


2 . 6 x 1 0 5 


1 . 86 xl 0 3 


12071 


8 . 0 % 


2.5x10 s 


1.79xl0 3 


12871* 


6.5% 


2 . 2 x 1 0 5 


1.52xl0 3 


135ti* 


5.5% 


1.9xl0 5 


1.33xl0 3 



Table 2. 300 A taper optimization results. 
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Table 2 summarizes the taper optimization results for a peak current of 300 A and 
an initial dimensionless optical amplitude of a 0 = 10. It was determined that for these 
parameters the optimum taper rate is 8 = 1 1271. This increased the efficiency and single- 
pass gain. The asterisk in the taper rate column indicates that the simulation included a 
hole was present in the center of the optical mode at the end of the undulator. 
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Figure 1 1 . 300 A peak current simulation with 8 = 1 12n. 

Figure 1 1 shows the results of the 3D diffraction simulation for the optimum taper 
rate determined from Table 2. Reducing the taper to 8 = 1 1271 allowed for continued 
guiding of the optical mode after the taper starts at x = 0.5. This is evident in the upper- 
left window which shows the majority of the optical power remains centered around the 
electron beam from x = 0 to x = 1. The upper-middle window also demonstrates improved 
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guiding, because no hole develops in the center of the optical mode. The bottom-middle 
window shows the bunched electrons trapped in the center of the phase-space plot. Note 
that around 50 % of the electrons remained trapped near resonance, v= 0, with the 
reduced taper of 8 = 1 1271. When the taper rate is too large, electrons cannot be captured 
by the optical field and are detrapped. 

The same optimization was performed for a peak current of 370 A realizing an 
optimal taper rate of 8 = 135 ti, with results summarized in Table 3. This provided only a 
slight increase in efficiency and single-pass gain over the design taper of 5 = 1 28 ti . The 
higher beam current allows for increased optical guiding which enables larger taper values 
to be used effectively. 

As can be seen in Figure 10, the optical mode has a hole in the final optical mode 
due to a loss of guiding of the optical mode along the undulator. Since the electron beam 
has a small radius of 0. 1 8mm, a large current density is necessary for guiding. While 



guiding occurs initially, when the taper turns on halfway through, the undulator diffraction 



8 




P(z) 


G f 


90ti 


8.2% 
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12071 


9.9% 
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130ti 
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4.0x1 0 5 


2.85xl0 3 


135ti 


10.3% 


4.1xl0 5 


2.86xl0 3 


140ti 


10.3% 


4.0xl0 5 


2.86xl0 3 


145ti 


10.0% 


4.0xl0 5 


2.86xl0 3 


150ti 


9.4% 


3.7xl0 5 


2.64x1 0 3 


155ti* 


8.2% 


3.3xl0 5 


2.37xl0 3 



Table 3. 370 A taper optimization results. 
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effects begin to take over. The optical mode begins to grow and the slope in the optical 
phase is negative, indicating a loss of guiding. In order to prevent this, the taper rate must 
not be too large for a given peak current. In the case of the 300 A peak current, a 
decreased taper rate did not result in a hole in the center of the optical mode leading to a 
corresponding increase in efficiency and power output. The change in taper rate allowed 
for increased optical guiding in the second half of the undulator. Higher currents allow the 
use of a higher taper rate before diffraction effects overcome optical guiding. This is 
evident when comparing Table 2 and Table 3. The hole in the optical mode does not 
occur until 5 = 1 55 tc for a peak current of 370 A, while at a lower current of 300 A the 
hole forms at a taper rate of only 8 = 128 tc. 

C. EFFECTS OF FEEDBACK 

Since this design is a cross between the oscillator and the amplifier, the amount of 
feedback could be essential in the understanding of the design. The RAFEL design uses 
only a small amount of feedback to seed the undulator. However, the amount of feedback 
for optimum efficiency has not been explored. In order to explore the effect of this aspect 
of the RAFEL, simulations varied the initial dimensionless optical field, ao, from 1 to 130. 
The single-pass 3D simulation is used as a representation of steady-state in an oscillator 
with the initial optical field determined by the amount of feedback. The inverse of the 
single-pass gain Gf is the amount of the feedback for high gain, or small feedback, 

/(G) = 1 / (l + G f ) . Therefore, the variation in feedback during the simulations ranged 

from about / = 1 0 _1 to / = 10” 5 . 
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The amount of feedback throughout this range surprisingly had almost no effect on 
efficiency and final power, as shown in Tables 4 and 5. The efficiency for the 300 A case 
was r| =7% ± 0.5%, while the efficiency for the 370 A case was again within 0.5% of the 
average value of r| = 10%. In both instances, saturation was reached because the output 



power was essentially the same due to the high gain at the beginning of the RAFEL 
undulator. 



ao 


*1 


P(x) 


G/ 


1 


7.2% 


2.2xl0 5 


1.6xl0 5 


5 


6.6% 


2.2x1 0 5 


6.1xl0 3 


10 


6.5% 


2.2x1 0 5 


1.5xl0 3 


40 


7.2% 


2.4x1 0 5 


106 


70 


6.9% 


2.4x1 0 5 


34 


100 


7.1% 


2.5x1 0 5 


17 


130 


7.5% 


2.7x1 0 5 


10 



Table 4. Feedback optimization for a 300 A peak current. 



ao 




p( t) 


G f 


1 


10.5% 


4.0xl0 5 


2.8xl0 5 


5 


9.9% 


3.9xl0 5 


l.lxlO 4 


10 


10.1% 


4.0xl0 5 


2.8xl0 3 


40 


10.7% 


4.3xl0 5 


191 


70 


10.4% 


4.3xl0 5 


61 


100 


10.2% 


4.3xl0 5 


30 


130 


10.2% 


4.4x1 0 5 


17 



Table 5. Feedback optimization for a 370 A peak current. 

Figures 12 and 13 correspond to a peak current of 300 A and initial dimensionless 
optical amplitude of a 0 = l and ao=130, respectively. These values cover the range of 
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feedbacks listed in Table 4, and represent feedback amounts on the order of / = 10 5 for 
Figure 12 and / = 10 -2 for Figure 13. In Figure 12, diffraction effects dominate over the 
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Figure 12. 3D simulation for 300 A and 10' 5 feedback, 
second half of the undulator resulting in a hole in the center of the optical mode. Also of 
note are the ineflfecient trapping of electrons, and the significant change in phase over the 
last half of the undulator. 

Figure 13 still shows the effects of diffraction of the optical mode along the 
undulator, but the majority of the mode is still centered around the electron beam at the 
end of the undulator with no hole developing in the optical mode. Notice that the phase 
change over the last half of the undulator is not as rapid, resulting in more efficient 
electron bunching. 
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Figure 13. 3D simulation for 300 A and 10' 2 feedback. 

When looking at the dimensionless power and efficiency in Tables 4 and 5, it 
appears that feedback does not play a significant role in the design. However, on closer 
examination it is evident that the shape of the optical mode is significantly effected. As the 
amount of feedback increases the hole in the optical mode diminishes and even goes away. 
This is important because the RAFEL laser output is the amount of optical power which 
exits the undulator from a hole in the center of the downstream mirror. The optical power 
which does not exit is fed back to seed the undulator. This means that only the outer 
portion of the optical mode is used for feedback. Since the hole in the output mirror has a 
diameter of 14 mm, only a tiny fraction is used for feedback. In fact, if the optical mode 
did not diffract sufficiently there would be no feedback at all. If the hole in the mirror 
were smaller a larger amount of the optical mode could be fed back. 
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D. EFFECTS OF INITIAL PHASE VELOCITY 



Because the RAFEL is actually an oscillator, the feedback used will allow the 
initial phase velocity to seek a maximum in gain. Since the initial phase velocity at weak 
optical fields is at v 0 = 0, the optical wavelength of the output could be noticibly different 
if the peak gain in strong fields is off resonance. In order to explore this effect numerous 
simulations like those in the previous sections were performed varying the optical 
amplitude a 0 and the initial phase velocity v c . 




Figure 14. Efficiency versus initial optical amplitude and phase velocity for 8 = 128n . 
Figure 14 shows the results of these simulations for a peak current of 300 A and a 
taper parameter 8 = 1 28rc . Constant in all these simulations was the dimensionless current 
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density j = 5.5 x 10 4 , the taper parameter 8 = 12871, and random Gaussian spread g g =13. 
The vertical axis is the efficiency, t|, and the other axes are the initial optical amplitude, a 0 , 
and the initial phase velocity, v 0 . As the optical amplitude increases, the efficiency peak 
drifts farther away from resonance culminating in a maximum efficiency of r| =9.3 % near 
a 0 = 20 and v 0 = 30, far off resonance. Over multiple passes ao will increase and v 0 will 
seek the maximum gain, the optical wavelength will change. The expected wavelength 
shift can be calculated from AX / X = Av / 2nN = 30 / 27tl 00 « 5% . Over the range 
simulated in Figure 14, the wavelength could shift by 5 %. This effect could be significant 
when trying to maximize power delivered through a small atmospheric window. 

E. EFFECTS OF PULSE EVOLUTION 

When short electron pulses are used they create short optical pulses. Over 
multiple passes, the arrival of these pulses at the beginning of the undulator must be 
synchronized. A single-pass simulation is useful in determining the possible short-pulse 
effects of a design. In order to run the simulations involving pulse evolution, an 
appropriate fillin g factor had to be chosen to describe the coupling in the transverse 
dimensions. Recall that the “filling factor” is the ratio of the area of the electron beam 
area to the optical beam area, F = r b / wl. The radius of the electron beam is known to be 
r b « 0.0 18 cm. The radius of the optical mode is less well defined due to the significant 
diffraction of the optical mode over the last meter of the undulator. Using the 3D 
simulations involved in the previous sections, it was evident that the optical mode over the 
untapered section of the undulator was relatively well behaved and could be estimated to 
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be about w 0 *0.027 cm. When the optical mode entered the tapered section, a significant 
portion of the mode diffracted away from the beam. For simplification a mode radius of 
w 0 *0.1 07 cm was chosen in the tapered section. This radius included that part of the 
optical mode having an optical amplitude of at least half of the maximum optical 
amplitude. Averaging these radii leads to an estimated filling factor of F * 0.07. 

Figure 15 below shows the results of a single-pass through the undulator. The 
dimensionless current density,/, in this simulation is actually jxF resulting in the value 
j = 4000. The initial phase velocity is v 0 = 0. The initial optical amplitude was a 0 = 10 
and the taper parameter matched the RAFEL design at 5 = 128n starting at x s = 0.5. The 
electron pulse length is given by a z = 3 . The upper-left window shows the optical pulse 
evolution along the undulator. The vertical axis for all three windows in the middle is 
dimensionless time evolving from t = 0 to x = 1 along the undulator. The horizontal axis 
represents the calculation window of width W = 6, which is larger than the electron pulse 
length plus slippage such that the optical pulse fits inside the window. The lower-left 
window shows the position of the short electron pulse. The light grey color represents the 
electron pulse at the beginning of the undulator, t = 0, and the black pulse represents the 
final position of the pulse after slippage at the end of the undulator, t = 1 . The upper- 
middle window shows the optical power spectrum development along the undulator. The 
development of small sidebands in this power spectrum is indicative of the trapped-particle 
instability. The upper-right window shows the electron phase velocity changes away from 
v = 0 along the undulator. The more electrons that are trapped, the larger the peak at 
v= 0. The amount the untrapped electrons drift out depends on the taper parameter taper 
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parameter 8 and the field ao. In this case, the taper starts halfway through the undulator 
so the untrapped electrons drift out to about 8 / 2 or phase velocity of v« 200. The 
bottom-middle window is a plot of the weak-field gain spectrum for reference. The lower- 
right window shows the optical power development along the undulator. The final gain 
for this simulation is G f = 4800 and the efficiency is t| = 4.9%. 




Figure 15. Single-pass pulse evolution simulation. 

The parameters in Figure 15 do not result in optimum efficiency. This can be seen 
by comparing the peaks in the upper-right window. Since about 50 % of the electrons 
should be trapped for optimum efficiency, the peaks in this window should have roughly 
equal areas. Too many electrons are untrapped; too few are traped near resonance v= 0. 
It is significant to point out that some indication of the trapped-particle instability is 
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evident even in this single-pass simulation. In order to verily the trapped-particle 
instability, multiple-pass simulations were performed using the same parameters. 

Short-pulse simulations were run over multiple passes using various values of 
desynchronization. As described previously, the desynchronization value, d, represents the 
timing of the optical feedback pulses relative to the incoming electron pulses. Figures 16, 
17, and 18 show the results of the multiple-pass pulse evolution simulations for 
desynchronizations, d = 0.1, d = 0.15, and d= 03 respectively. The feedback in these 

simulations is / = 1 / (l + G) where gain is G = {oj / a\ ) - 1, so that the feedback is 
/ = Oq / a 2 f . In these simulations, a 0 = a f exp(-l / 2Q) so / = exp(-l / Q). The smallest 
value of feedback at which the FEL would effectively operate was found to be Q = 0.15, 
corresponding to feedback / = 10~ 3 . The simulations would not reach saturation at lower 
feedback values because the taper rate gives low gain in weak fields reducing coupling 
between the electrons and the optical wave. 

The upper-left window represents the optical pulse evolution over n passes where 
black and dark grey indicates higher optical amplitude. The upper-middle window 
represents the power spectrum evolution over n-\ 00 passes. The upper-right window 
shows the electron distribution over n passes. The lower-left and lower-middle windows 
are the same as before, while the lower-right window now represents the total pulse power 
evolution over n passes. 

Of particular note in these multiple-pass simulations is the power spectrum in the 
upper-middle window of the simulations. Close examination shows that the peak in the 
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power drifts off resonance away from v = 0. This drift was anticipated in earlier 
simulations in the previous section. The simulations in this section correspond to an initial 
optical amplitude at the upper limit of Figure 14. The final optical amplitude is 
approximately a f =700 leading to a 0 = a f exp(-l / 2 Q) ~ 25 for Q = 0.1 5 . Measurement 

of the phase velocity deviation from resonance in Figure 16 showed a change in phase 
velocity of Av« 25, and a corresponding wavelength variation given by 
AX. / X = Av / 2nN = 25/ 2n \ 00 « 4%. This is in agreement with predictions in the 
previous section where peaks in the gain spectrum occur off-resonance by Av « 25. The 
optical power has followed the peaks in the gain spectrum in strong optical fields. 



**** FEL Pulse Evolution **** 
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Figure 16. Multiple-pass pulse evolution ( d = 0.1) simulation. 
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In Figure 16, the optical pulse is shortener with smaller values of 
desynchronization [3]. Additionally, the optical pulse is distorted showing the spikes 
indicating the presence of trapped-particle instability. The power spectrum broadens over 
multiple passes and the output power fluctuates significantly even after saturation is 
achieved. The power actually drops off after 30-40 passes and then reaches steady state. 
The final power at d = 0.1 is about half the power for optimum desynchronization. 



**** FEL Pulse Evolution **** 
j =4000 c z =3 Q=0. 15 c G =l . 3 d=0.15 
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Figure 17. Multiple-pass pulse evolution ( d = 0.15) simulation. 

Figure 17 also shows spikes in the optical pulse and broadening of the power 
spectrum which indicate the trapped-particle instability. At this intermediate value of 
desynchronization, the optical pulse is much wider and some electrons still remain trapped. 
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Figure 1 8 shows the effects of the trapped- particle instability with a broader power 
spectrum. 



**** FEL Pulse Evolution **** 
j =4000 c z =3 Q=0.15 (jg=1.3 d=0 . 3 
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Figure 18. Multiple-pass pulse evolution ( d = 03) simulation. 

Figure 19 shows a desynchronization plot for the RAFEL. The plot was compiled 
using a series of muliple pass pulse evolution simulations using the average power after 
saturation was reached. The average power was used because power fluctuated by as 
much as 20% after saturation was reached. The same parameters were used for these 
results as those for the previous multiple-pass pulse evolution simulations. The x’s 
indicate simulation results and the line traces the estimated curve shape. The o' s 
correspond to the simulations shown in Figures 16, 17, and 18. The trapped-particle 
instability caused the fluctuations in power over the entire desynchronization range. The 
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high gain of the RAFEL increases the values of d beyond those associated with traditional 
FEL oscillators, but the curve shape remains the same. 




Figure 19. Desynchronization plot for the RAFEL. 



F. IDEAL RAFEL OPERATION 

In order to get a more complete understanding of the operation of the RAFEL 
design, one dimensional simulations were run with no diffraction and a perfect electron 
beam. Figure 20 shows the expected small signal gain spectrum using a 0 = 1 x 10 -6 with 
j = 4000 and taper rate 8 = 12871. The upper window is the gain spectrum for initial 
phase velocities, -60 < v 0 < 30. The peak of the gain curve occurs at v 0 « 0 with a 
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maximum of G =41,280. The lower window is a plot of the optical phase versus initial 
phase velocity. 




Figure 20. Small signal gain spectrum. 




**** Gain and Phase Curves **** 
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Figure 21. Gain spectrum for a 0 = 30. 
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Figure 21 shows the gain spectrum for a 0 = 30 with the same beam and taper 
parameters as the previous figure. The gain spectrum has a distinctive shape and 
corresponds to a situation in which a larger amount of optical power has been fed back to 
the beginning of the undulator. The main peak has shifted to a maximum at v 0 = 203 . 
Recall that Figure 14 showed a peak in the efficiency around this initial phase velocity for 
a 0 - 25. The more dramatic change in the gain spectrum is the development of significant 
gain peaks for v 0 < -30. This can be explained using phase space plots. 




Figure 22. Small signal phase space plot for LANL design. 

Phase space plots are useful in understanding the physics behind the gain curve 
results. Specifically where the electrons are bunching and how many are getting trapped 
determines the resulting gain. Figure 22 is a phase space plot at v 0 = 0 using the same 
parameters as Figure 20. The main window shows the phase space plot with electron 
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phase velocity along the vertical axis and the electron phase along the horizontal axis. It is 
evident in the plot that none of the electrons were trapped. This is not surprising since the 
optical amplitude at = 0.5 is a = a 0 -fG » 10” 4 and is not large enough to take advantage 
of the tapered undulator. The upper-right window shows the natural log of the gain along 
the undulator. The gain grew significantly during the untapered first half of the undulator, 
but then leveled out over the last meter where the undulator was tapered but no electrons 
were trapped. This is indicative of the ineffectiveness of the taper at small initial optical 
amplitude when there is no trapping. The lower-right window shows the development of 
the optical phase along the undulator. The phase shift was about n over the first half of 
the undulator and fluctuated around this value over the remainder of the undulator. 



*** FEL Phase Space Frolatian *** 
j =4000 a o =30 v o =20.5 N=100 
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Figure 23. Phase space plot for a 0 = 30. 
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Figure 23 is a phase space simulation at the more optimum value of v 0 = 20.5 
corresponding to the peak in Figure 21 . This shows the electron phase development for 
this optimum condition. About 50% of the electrons are trapped near resonance v = 0 and 
bunched around C, = n. The taper works in this case because the field at x s = 0.5 is 

a = a 0 4 g « 10 3 , and is sufficiently strong to trap electrons. The gain continues to grow 
throughout the undulator length indicating the taper was effective for this larger optical 
amplitude. 

In Figure 20 the gain spectrum is relatively close to that of an untapered undulator 
because none of the electrons are trapped. Since the taper is ineffective, the gain spectrum 
resembles that of an untapered undulator with half the length of the actual undulator. The 
unique shape of the gain spectrum in Figure 21 is a result of the taper starting point. The 
peaks in the gain spectrum corresponding to negative initial phase velocities occur because 
as the electrons increase in phase velocity from the tapers artificial acceleration they 
become trapped in the separatrix when it is centered near C, + <|) = n allowing them to 
contribute to the optical mode. The dips near these peaks correspond to electrons being 
trapped when the separatrix is not in the optimum area of phase space for gain. The taper 
should be turned on when the electrons are giving up energy to the optical mode and give 
an artificial acceleration to keep them giving energy to the optical mode. If the taper is 
turned on when the electrons are taking energy from the optical mode, it is not as effective 
or efficient. Also, if the taper amount is too high or too low there will be a change in 
optical phase and the separatrix will drift across phase space causing the electrons to leak 
out and diminishing gain. 
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V. CONCLUSION 



The U.S. Navy needs to develop a “speed of light” missile defense capability as an 
alternative to current technology. The FEL has the potential to provide this capability. 

The LANL RAFEL design advances research into this capability by combining the 
oscillator and amplifier FEL characteristics. In order to better understand the 
characteristics of the RAFEL design, simulations were done addressing feedback and taper 
optimization, initial phase velocity effects, short-pulse effects, and ideal RAFEL operation. 

The taper optimization was performed for peak currents of 300 A and 370 A and a 
taper starting point halfway through the undulator. In the 300 A case, peak efficiency can 
be achieved by decreasing the taper rate from the current design of 6 = 128rc down to 
5 = 1 12n or 26% taper instead of 30% taper. For the higher peak current of 370 A, the 
optimum efficiency occurred at 5 = 135 tc. The key to these results is that the higher 
current led to more effective optical guiding. Higher taper rates cause optical diffraction 
to overcome optical guiding causing a hole to develop in the center of the optical mode. 

The feedback optimization was performed using the designed taper rate of 
8 = 128n for the same peak currents. This was accomplished by varying the initial optical 
amplitude corresponding to feedbacks of / = 10‘ 5 to 10‘ 2 . In the 300 A and 370 A cases, 
the efficiency was not significantly effected over this range. However, when the initial 
optical amplitude was reduced further, a 0 < 0.1, there was no hole in the center and the 
efficiency dropped dramatically. Additionally, over the entire range of feedbacks the shape 
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of the optical mode varied. For small feedback there was a significant hole in the center of 
the optical mode, while at larger feedback there was no hole. Since the feedback amount 
is dependent on the shape of the optical mode, this introduces another variable in the 
RAFEL design. Until a steady-state mode of operation is reached, the amount of feedback 
will vary. A smaller aperture in the downstream mirror would increase feedback and 
possibly increase efficiency during start up. 

The phase velocity giving peak efficiency depends significantly on the initial optical 
amplitude. Since the RAFEL design involves feedback, the dependence of efficiency on 
the initial phase velocity could cause a change in the optical wavelength. The wavelength 
change would be on the order of a few percent, and can be taken into account when 
designing a RAFEL to operate within a desired wavelength window. 

Since the electron beam is really a series of short electron pulses, the effects of 
short pulses were addressed in simulations. For multiple-pass short -pulse simulations, the 
RAFEL design would only operate effectively at feedbacks greater than / = 10” 3 . These 
simulations showed that the optical pulse was significantly distorted by the trapped- 
particle instability, and the output power fluctuated after reaching saturation. 

Simulations indicated that the RAFEL design depends significantly on the taper 
starting point and taper rate. Since the RAFEL starts from weak fields and tapering is 
only effective in strong optical fields, it is inherently inefficient at start up. While the gain 
over the first meter of the undulator is high, it is not high enough to allow efficient use of 
tapering when the initial optical amplitude is too small. The feedback loop enables the 
RAFEL design to get through this inefficient phase by seeding the undulator with a small 
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amount of optical power. The result is a gain spectrum that varies significantly with the 
initial optical amplitude. As the amount of feedback increases, the fixed taper rate varies 
in effectiveness. Like the 3D simulations involving initial phase velocity, the ideal 
simulations showed a shift in the gain peak toward larger values of v 0 as the initial optical 
amplitude was increased. 

The RAFEL design effectively combines the attributes of the oscillator and 
amplifier FEL designs. It is inherently inefficient at the beginning of each macropulse 
when optical amplitudes are small. However, saturation and efficient operation can still be 
reached due to the combination of high gain and feedback. Detailed experiments, or 
simulations, should be done on the feedback system to better understand the RAFEL 
operation, and optical mode development. One of the advantages of the RAFEL design is 
reduced power on the mirrors. As future designs work toward higher power, better 
understanding the optical mode development over multiple passes will be important to 
maintain this advantage. Future work could involve detailed simulation of the feedback 
loop or design modifications. 
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